MEASUREMENT OF -> £+l> AND THE DECAY CONSTANT F Ds 
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I report preliminary CLEO-c results on purely leptonic decays of the D s using 195 pb 1 of data 
at 4.170 GeV. We measure fo 3 = 280.1 ± 11.6 ± 6.0 MeV, and / + /f D+ = 1.26 ± 0.11 ± 0.03. 
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1. Introduction 

To extract precise information from B — B 
mixing measurements the ratio of "leptonic 
decay constants," /j for Bd and B s mesons 
must be well known. 1 Indeed, the recent mea- 
surement of B° — B s mixing by CDF 2 has 
pointed out the urgent need for precise num- 
bers. The fi have been calculated theoreti- 
cally. The most promising of these calcula- 
tions are based on lattice-gauge theory that 
include the light quark loops. 3 In order to 
ensure that these theories can adequately 
predict fB s /fB d it is critical to check the 
analogous ratio from charm decays f D + / fo+ ■ 
We have previously measured fjj+ . 4 ' 5 Here 
I present the most precise measurements to 
date of f D + and f D +/f D+ . 

In the Standard" Model (SM) the D s me- 
son decays purely leptonically, via annihi- 
lation through a virtual W + , as shown in 
Fig. 1. The decay width is given by 6 



Fig. 1. The decay diagram for — > i + v. 

would simply change the value of /j extracted 
using Eq. (1). We can, however, measure the 
ratio of decay rates to different leptons, and 
the predictions then are fixed only by well- 
known masses. For example, for t + v to [i + i>: 
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where m g + and M D + are the £ + and D + 
masses, \V CS \ is a CKM element equal to 
0.9737, and Gf is the Fermi constant. 

New physics can affect the expected 
widths; any undiscovered charged bosons 
would interfere with the SM W + . These ef- 
fects may be difficult to ascertain, since they 



T(D -> n+v) 



(2) 

Any deviation from this formula would 
be a manifestation of physics beyond the SM. 
This could occur if any other charged inter- 
mediate boson existed that coupled to lep- 
tons differently than mass-squared. Then the 
couplings would be different for muons and 
t's. This would be a clear violation of lepton 
universality. 7 

2. Experimental Method 

In this study we use 195 pb -1 of data pro- 
duced in e + e~ collisions using the Cornell 
Electron Storage Ring (CESR) and recorded 
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near 4.170 GeV. Here the cross-section for 
D*+D~+D+D*- is -1 nb. We fully recon- 
struct one D s as a "tag," and examine the 
properties of the other. In this paper we des- 
ignate the tag as a Dj and examine the lep- 
tonic decays of the , though in reality we 
use both charges. Track selection, particle 
identification, 7, 7r°, K$ and muon selection 
cuts are identical to those described in Ar- 
tuso et al. 

The D~ decay modes used for tagging 
are listed in Table 1. The number of signal 
and background events are determined by fits 
to the invariant mass distributions. 

Table 1 . Tagging modes and numbers of signal and 
background events, within ±2.5cr for all modes, ex- 
cept rjp + (±2cr), from two-Gaussian fits to the in- 
variant mass plots. 



Mode 


Signal 


Background 


K+K--K- 


8446 ± 160 


6792 


K S K~ 


1852±62 


1021 


jpx~ ; r) — » 77 


1101 ±80 


2803 


J]'tT~; Vj — > 7T ^ 7T Tj 


786±37 


242 


4>P~ 


1140±59 


1515 


7r + 7r _ 7r — 


2214±156 


15668 


K*-K*° 


1197±81 


2955 


r)p- 


2449±185 


13043 


Sum 


19185 ± 325 


44039 



We search for three separate decay 
modes: Df —> /J, + is, Df — > t + v. t + — > tt + V 
and t + — > e + vV. For the first two analyses 
we require the detection of the 7 from the 
D* — > jD s decay. Regardless of whether or 
not the photon forms a D* with the tag, for 
real D*D S events, the missing mass squared 
recoiling against the photon and the D~ tag 
should peak at the mass and is given by 

MM* 2 = {Ecu -E D - E y f - (~p2 - , 

where Eqm is the center of mass energy, Ejy 
(pd) and E 1 (p^) are the energy (momen- 
tum) of the fully reconstructed D~ tag, and 
the additional photon. In performing this 
calculation we use a kinematic fit that con- 
strains the decay products of the D~ to the 



known D s mass and conserves overall mo- 
mentum and energy. 

The MM* 2 from the Dj tag sample data 
is shown in Fig. 2. Fitting shows a yield of 
12604±423 signal events. Restricting to the 
interval 3.978 >MM* 2 > 3.776 GeV 2 , we are 
left with 11880±399 events. The systematic 
error is ±4.3%. 
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Fig. 2. The MM* 2 distribution from events with 
a photon in addition to the tag. The curve is 
a fit to the Crystal Ball function and a 5th order 
Chebychev background function. 

Candidate Df — > fi + u events are 
searched for by selecting events with only 
a single extra track with opposite sign of 
charge to the tag; we also require that there 
not be an extra neutral energy cluster in ex- 
cess of 300 MeV. Since here we are searching 
for events where there is a single missing neu- 
trino, the missing mass squared, MM 2 , eval- 
uated by taking into account the seen 
D~ , and the 7 should peak at zero, and is 
given by 

MM 2 = (E CM — Ed — E 1 — E„f (3) 
- (-Pd - ih - p£) , 

where E^ (p^) is the energy (momentum) of 
the candidate muon track. 

We also make use of a set of kinematical 
constraints and fit the MM 2 for each 7 can- 
didate to two hypotheses one of which is that 
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the Dj tag is the daughter of a D*~ and the 
other that the D* + decays into -fDf, with 
the subsequently decaying into [i + v. The 
hypothesis with the lowest x 2 is kept. If there 
is more than one 7 candidate in an event we 
choose only the lowest \ 2 choice among all 
the candidates and hypotheses. 

The kinematical constraints are the to- 
tal momentum and energy, the energy of the 
either the D* or the D s , the appropriate 
D* — D s mass difference and the invariant 
mass of the D s tag decay products. This 
gives us a total of 7 constraints. The missing 
neutrino four-vector needs to be determined, 
so we are left with a three-constraint fit. We 
preform a standard iterative fit minimizing 
X 2 . As we do not want to be subject to sys- 
tematic uncertainties that depend on under- 
standing the absolute scale of the errors, we 
do not make a \ 2 cu t, but simply choose the 
photon and the decay sequence in each event 
with the minimum \ 2 ■ 

We consider three separate cases: (i) the 
track deposits < 300 MeV in the calorime- 
ter, characteristic of a non-interacting 7r + or 
a (ii) the track deposits > 300 MeV in 
the calorimeter, characteristic of an interact- 
ing 7r + ; (hi) the track satisfies our e + selec- 
tion criteria. Then we separately study the 
MM 2 distributions for these three cases. The 
separation between \i + and 7r + is not unique. 
Case (i) contains 99% of the fi + but also 60% 
of the 7r + , while case (ii) includes 1% of the 
fj+ and 40% of the tt+. 5 

The overall signal region we consider is 
below MM 2 of 0.20 GeV 2 . Otherwise we ad- 
mit background from 7]ir + and K°tt + final 
states. There is a clear peak in Fig. 3(i), due 
to D+ — > n + v. Furthermore, the events in 
the region between ijl^v peak and 0.20 GeV 2 
are dominantly due to the t + v decay. 

The specific signal regions are defined as 
follows: for fj+u, 0.05 >MM 2 > -0.05 GeV 2 , 
corresponding to ±2tr or 95% of the signal; 
for tv, t+ -> 7T+77, in case (i) 0.20 >MM 2 > 
0.05 GeV 2 and in case (ii) 0.20 >MM 2 > 
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Fig. 3. The MM 2 distributions from data using DJ 
tags and one additional opposite-sign charged track 
and no extra energetic showers (see text), (a) Case 
(i) the single track deposits < 300 MeV of energy in 
the calorimeter. The peak near zero is from — t 
fj, + u events, (b) Case (ii) Track deposits > 300 MeV 
in crystal calorimeter but is not consistent with being 
an e~K (c) Case (iii) the track is identified as an e+. 

—0.05 GeV 2 . In these regions we find 64, 24 
and 12 events, respectively. The correspond- 
ing backgrounds are estimated as 1, 2.5 and 1 
event, respectively. The branching fractions 
are summarized in Table 2. The absence of 
any detected e + opposite to our tags allows 
us to set the upper limit listed in Table 2. 

The — > r + f, r + — » e + vv mode is 
measured by detecting electrons of opposite 
sign to the tag in events without any addi- 
tional charged tracks, and determining the 
unmatched energy in the crystal calorimeter 
(Eg^ ra ). This energy distribution is shown 
in Fig. 4. Requiring Egg ra < 400 MeV, en- 
hances the signal. The branching ratio re- 
sulting from this analysis is listed in Table 2. 
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Table 2. Measured DJ" Branching Fractions 



Final State 



B (%) 









0.657 ± 0.090 ± 0.028 








0.664 ±0.076 ±0.028 


■H 


-i/, (r+ - 


+ 7T+1/) 


7.1 ±1.4 ±0.3 




"I/, (r+ - 




6.29 ±0.78 ±0.52 


-H 


~v (averaj 


;e) 


6.5 ±0.8 


e 4 


V 




< 3.1 x 10" 4 (90% cl) 



f From summing the fi^u and t + u contribu- 
tions for MM 2 < 0.20 GeV 2 . 



I ' ' ' ' I ' ' ' ' 

° Data 

— MC Total 

— MC Signal 
MC Background 

MC K°eV 




E^ ra (GeV) 

Fig. 4. The extra calorimeter energy from data 
(points), compared with the Monte Carlo simulated 
estimates of semileptonic decays in general (dotted), 
the K°e^"u mode specifically (shaded), as a sub-set 
of the semileptonics, and the expectation from sig- 
nal (dashed). The peak near 150 MeV is due to the 
7 from D* — > *yD s decay. (The sum is also shown 
(line).) The arrow indicates the selected signal region 
below 0.4 GeV. 

3. Conclusions 

Lepton universality in the SM requires that 
the ratio R from Eq. 2 be equal to a value of 
9.72. We measure 

T(D+ - t+u) 



R 



T(D+ -» fi+v) 



9.9 ± 1.9 . 



(4) 



Thus we see no deviation from the predicted 
value. Current results on D + leptonic de- 
cays also show no deviations. 8 Combining 
all three branching ratios determinations and 



f+ 

J D 



222.6 ± 16.7±|;| MeV, 



using t d + =0.49 ps to find the leptonic width, 
we find 

f D , = 280.1 ±11.6 ±6.0 MeV. (5) 
Using our previous result 

(6) 

provides a determination of 

f D +/fu+ =1-26 ±0.11 ±0.03. (7) 

These preliminary results are consis- 
tent with most recent theoretical models. 
As examples, unquenched lattice 9 predicts 

I. 24 ±0.01 ±0.07, while one quenched lattice 
calculation 10 gives 1.13 ± 0.03 ± 0.05, with 
other groups having similar predictions. 11 

Acknowledgments 

This work was supported by the National 
Science Foundation. I thank Nabil Menaa 
for essential discussions. 

References 

1. G. Buchalla, A. J. Buras and M. E. Lauten- 
bacher, Rev. Mod. Phys. 68, 1125 (1996). 

2. A. Abulencia et al. (CDF), "Observation 
of B S B S Oscillations," [hep-ex/0609040]; 
see also V. Abazov et al. (DO), [hep- 
ex/0603029]. 

3. C. Davies et al, Phys. Rev. Lett. 92, 022001 

(2004) . 

4. M. Artuso et al. (CLEO), Phys. Rev. Lett. 
95, 251801 (2005) [hep-ex/0508057]. 

5. G. Bonvicini, et al. (CLEO) Phys. Rev. 
D70, 112004 (2004) [hep-ex/0411050]. 

6. D. Silverman and H. Yao, Phys. Rev. D38, 
214 (1988). 

7. J. Hewett, "Seaching For New Physics with 
Charm," SLAC-PUB-95-6821 (2005) [hep- 
ph/9505246]; W.-S. Hou, Phys. Rev. D48, 
2342 (1993). 

8. P. Rubin et al. (CLEO), Phys. Rev. D73, 
112005 (2006). 

9. C. Aubin et al, Phys. Rev. Lett. 95, 122002 

(2005) . 

10. T. W. Chiu et al, Phys. Lett. B624, 31 
(2005) [hep- ph/0506266]. 

II. See references to other theoretical predic- 
tions in [4]. 



